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Abstract 

Standard Schramm-Loewner evolution (SLE) is driven by a con- 
tinuous Brownian motion which then produces a trace, a continuous 
fractal curve connecting the singular points of the motion. If jumps 
are added to the driving function, the trace branches. In a recent pub- 
lication [1] wc introduced a generalized SLE driven by a superposition 
of a Brownian motion and a fractal set of jumps (technically a stable 
Levy process). Wc then discussed the small-scale properties of the re- 
sulting Levy-SLE growth process. Here we discuss the same model, 
but focus on the global scaling behavior which ensues as time goes to 
infinity. This limiting behavior is independent of the Brownian forc- 
ing and depends upon only a single parameter, a, which defines the 
shape of the stable Levy distribution. We learn about this behavior 
by studying a Fokker- Planck equation which gives the probability dis- 
tribution for endpoints of the trace as a function of time. As in the 
short-time case previously studied, wc observe that the properties of 
this growth process change qualitatively and singularly at a = 1. We 
show both analytically and numerically that the growth continues in- 
definitely in the vertical direction for a > 1, goes as \ogt for a ~ 1, and 
saturates for a < 1. The probability density has two different scales 
corresponding to directions along and perpendicular to the boundary. 
In the former case, the characteristic scale is X{t) ~ t^^". In the lat- 
ter case the scale is Y{t) ^ A + Bt^-^^°' for a ^ 1, and Y(t) ~ Int 
for a = I. Scaling functions for the probability density are given for 
various limiting cases. 
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1 Introduction 



The study of random conformally-invariant clusters that appear at critical 
points in two-dimensional statistical mechanics models has been made rig- 
orous with the invention of the so-called Schramm-Loewner evolution (SLE) 
[^. SLE refers to a continuous family of evolving conformal maps that spec- 
ify the shape of a part of a critical cluster boundary. By now SLE has 
been justly recognized as a major breakthrough, and there are several re- 
view papers and one monograph devoted to this beautiful subject, see Refs. 

SLE describes a curve, called trace, growing with time from a boundary 
in a two-dimensional domain which is usually chosen to be the upper half 
plane. SLE is based on the Loewner equation in which the shape of the 
growing curve is determined by a function of time ^(t) which in SLE is 
taken to be a scaled Brownian motion. Such a choice of the driving function 
produces continuous stochastic, fractal and conformally invariant curves, 
— the kind that appears as the scaling limit of various interfaces in many 
two-dimensional critical lattice models and growth processes of statistical 
physics. Well-known examples include boundaries of the Fortuin-Kastelyn 
clusters in the critical g-state Potts model, loops in the 0(n) model, self- 
avoiding and loop-erased random walks. 

In Ref. Q we generalized SLE to a broader class for which ^(t) is a 
Markov process with discontinuities. More specifically, we have studied the 
Loewner evolution driven by a linear combination of a scaled Brownian mo- 
tion and a symmetric stable Levy process. The growing curve then exhibits 
branching. This generalized process might be useful to describe many tree- 
like growth processes, such as branching polymers and various branching 
growth processes which evolve in time. 

Such generalized SLEs driven by Levy processes (Levy-SLE for short) 
have also been of interest to mathematics community. Our results [Q] on 
various phase transitions in Levy-SLE have been put on rigorous basis in Ref. 



and further properties have been studied in Refs. [|T^, 13]. The interest 
of mathematicians in these Levy-SLE processes is partially motivated by the 
sugg estion ||, |l|] that they may produce fractal objects with large values 
of multifractal exponents for harmonic measure. Harmonic measure can 
be thought of as the charge distribution on the boundary of a conducting 
cluster. On fractal boundaries such a distribution is a multifractal, and 
in the case of critical clusters (whose boundaries are SLE curves) the full 
spectrum of multifractal exponents has been obtained analytically, see Refs. 



4|, 15, |l^, |r^, 18 1 for various derivations and discussion. 
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While our previous paper ||] focused on local properties of Levy-SLE, 
here we study the global behavior of the growth in the upper half plain. The 
present paper is structured as follows. In Section |2| we define our model and 
briefly state our previous results on phase transitions in the local behavior 
of the model. We also present our new results on the global behavior of 
Levy-SLE. In Section |3| we derive the Fokker-Planck equation governing the 
evolution of the probability distribution for the tip of the Levy-SLE. The 
equation is our main tool for analysis of the long time global behavior of 
the growth. We give a qualitative description of the growth and explain the 
approximations that go into the solution of the Fokker-Planck equation in 
Section ^. Actual solution of the Fokker-Planck equation and comparison 
with results from numerically calculated trajectories is given in Section ^. 
We conclude in Section |6|. Some technical details are presented in Appen- 
dices. 



2 The model and the results, old and new 

Loewner evolution is a family of conformal maps that appears as the solution 
of the Loewner differential equation (see, for example, Ref. ^ for details) 

2 

dtgt{z) = —— — , go{z) = z. (1) 

9t[z)-C{t) 

valid at any point z in the upper half plane until (and if) this point becomes 
singular at some (possibly infinite) time r^: ^(t^) = gr^{z). The set of all 
singularities is called the hull and the point at which the hull grows is called 
the tip. The tip 7(t) is defined via its image ^{t) = gt{'y{t)). More formally, 

7(i)= lim 9^\w), (2) 

where the limit is taken in the upper half plane. The trace is the path left 
behind by the tip (the existence of the trace in the setting of this paper 
has been shown in Ref. |]l2|). The shape of the growing trace (and the 
hull) is completely determined by the driving function ^(t). At any time 
the function gt{z) conformally maps the exterior of the growing hull to the 
upper half plane, see Fig. [l|. We refer to the z plane where the growth 
occurs as "the physical plane", and to the w plane as "the mathematical 
plane" . 

Naturally, if £,{t) is a stochastic process the shape of the growing trace is 
also stochastic. The growth process is then a stochastic (Schramm-) Loewner 
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Figure 1: The Loewner evolution shown for the case when the growing hull 
is a smooth curve. The complement of the segment of the curve (up to its 
tip 7(t)) in the "physical" z plane is mapped to the entire upper half of the 
"mathematical" w plane by the function gt{z). 



evolution (SLE). The standard SLE has a driving function ^{t) = y/K.B{t), 
where B{t) is a normalized Brownian motion and k > is the diffusion 
constant. Many important properties of this process have been established 
inRef. [||]. 

In Ref. mi we have generalized SLE to 

e(t) = 0^i?(t) + ci/"L„(t), (3) 



where La{t) is a normalized symmetric a-stable Levy process |2^, 21, p^ , 
and c > is the "diffusion constant" associated with it. The process La{t) 
is composed of a succession of jumps of all sizes. Unlike a Brownian motion, 
La{t) is discontinuous on all time-scales. Therefore, the addition of a Levy 
processes to the driving force of SLE introduces branching to the trace. 

The probability distribution function of c^^"La{t) is given by the Fourier 
transform 

p{x,t)= r ^^-ikx^-ct\k\-^ 

As it is known in the theory of stable distributions [^] , only for < a ^ 2 
this Fourier transform gives a non-negative probability density. For < a < 
2 the function P{x, t) decays at large distances as a power law: 

ct 

x^oo: Pix,t) ^ I (^) 
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so that the process scales as 



(|L,(t)|^)oct^/- (6) 

for any 5 < a. For 5 ^ a this average is infinite. For q = 2 the process 
L2{t) is the standard Brownian motion B{t) and P{x,t) is Gaussian. 

We studied the short-distance properties of the Levy-SLE process in Ref. 
[Ql . At short times and distances the process is dominated by the Brownian 
motion and the deterministic drift term (see Eq. ([T^)), whereas at long 
times it is dominated by Levy flights. The crossover between short and long 
time behavior happens at the time 

/l\l/(2-a) 

^0 - (^) . (7) 

This also defines a spatial crossover at length scales Iq oc ^/tQ. For scales 
smaller than Iq the trace behaves like standard SLE, while for scales much 
larger than Iq it spreads in the x direction forming tree-like structures. 

In our previous paper Q, using both analytic and numerical considera- 
tions, we determined the probability that a point on the x axis is swallowed 
by the trace. The trace shows a qualitative change in its small-distance, 
small-time behavior as k and a each pass though critical values, respec- 
tively at four and one. The transition at k = 4 is quite analogous to the 
known transition of standard SLE For the new transition at a = 1, the 
trace forms isolated trees when a < 1 or a dense forest when a > 1. 

The latter phase transition at a = 1 was recently studied rigorously in 



Ref. 1 11 1 which expanded the implications of the phase transition to the 
whole plane at the limit t ^ oo. For k > 4 a point in the upper half plane 
is swallowed almost surely for a > 1, while it is swallowed with probability 
smaller than one for a < 1. For k < 4 and < a < 2 the swallowed points 
on the plane form a set of measure zero. 

The large-scale implications of the a = 1 transition can be seen in Figure 
^ which shows the shape of the trace at long times. For a < 1 the stochastic 
evolution produces isolated tree-like structures which are limited in height. 
For a > 1 the evolution produces an "underbrush" in which structures pile 
on one another and thereby continue to increase their height. 

In the rest of the paper we establish the following. The growth at long 
times is characterized by two very different length scales X(t) and Y{t) (with 
X{t) ^ yit)) which can be thought of as the typical size of the growing hull 
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Figure 2: Examples of traces produced by Levy-SLE at long times, up to 
t = 9000. For the first, second and last thirds of the time interval, the traces 
are correspondingly colored red, blue, and yellow. Top panel: a = 0.7. The 
trace looks like many isolated trees whose height saturates at long times. 
Bottom panel: a = 1.5. Now the tip of the growing trace keeps landing on 
the previously grown "bushes" so that the trace extends indefinitely in the 
vertical direction as time increases. Notice the difference in scales for the 
y axis between the two panels, as well as the much larger spread in the x 
direction in the top panel. The trace was produced using stable Levy forcing 
and c = 10, time step r = 10~^. The trace was calculated only at times 
when the forcing makes a large jump dC > V2OO T. 
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in the x and y directions. More specifically, we find that 

X{t) ~ t^/", < a < 2, (8) 

+ (9) 

I Int, a = 1. 

(The constants A and B depend upon a.) These scales enter the scaling form 
of the joint probability distribution p{x, y, t) for the real and imaginary parts 
of the tip 7(t) of the Levy-SLE, for which we give explicit results in various 
limiting cases in Section ^, where we also compare analytical results with 
extensive numerical simulations. 



3 Derivation of the Fokker-Planck equation 

We are interested in characterizing the probability distribution for the point 
7(t) at the tip of the trace in the ensemble provided by different realizations 
of the SLE stochastic process. Eq. (|2|) implies then that we should study the 
inverse map g'^^ . However, this is rather difficult, since the map g^^ satisfies 
a partial differential equation instead of an ODE. There is a way out which 
is rather well known and has been successfully used before [13, 14, 19|. It 
happens that one needs to consider the backward time evolution: 

= - h{u,)-i{ty /oH = ^- (10) 

The relation of the original Loewner evolution (||) and the backward one 
( p!o| ) in the stochastic setting is as follows. If ^(t) is a symmetric (in time) 
process with independent identically distributed increments, which is the 
case for a Levy process, then it is easy to show that for any fixed time t the 
solution ft{w) of the backward equation ([lO|) has the same distribution as 

gt^{w — C{t)) + ^{t), see Refs. [13, O]. Using the symbol = for equality of 



distributions for random variables, we can write 

f,{w)^grHw-m)+m- (n) 

It is useful to introduce a shifted conformal map 

ht{z)=gt{z)-C{t), (12) 
for which the Loewner equation acquires the Langevin-like form: 

dthtiz) = - mt), ho{z) = z, (13) 
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assuming that ^ vanishes at t = 0. The first term is a deterministic drift 
and the second — a random noise. The tip j{t) is now mapped to zero, and 
this can be taken as the definition of the tip. More formally, 

7(t)= lim V'M (14) 

where the limit is taken in the upper half plane. 

In terms of the shifted map the equality of distributions (11) can be 
written as 

ft{w) - m = ht\w). (15) 

The left hand side Zt = ft{w) — i{t) of this equation satisfies the Langevin- 
like equation 

2 

dtzt = dti{t), zq = w, (16) 

and in particular, if we set w = in this equation, the resulting stochastic 
dynamics should be the same as that of the tip of the trace 7(t). 

Before we convert the Langevin-like equation ( p!6|) to our main analytical 
tool, the corresponding Fokker-Planck equation, let us review again the 
correspondence between the forward and backward flows and illustrate it 
with figures. Equation ( [T^ ) describes a flow in which wt = ht{z) follows a 
trajectory of a particle in the w plane, z being its initial position. Separating 
the real and imaginary parts of wt = ut + ivt, we get a system of coupled 
equations 

dtut = n^"* n - dtC{t), uo = X, 

uf + vf 

^t'"i = 2~7^' ^0 = y, (17) 

describing such a trajectory. As in many modern versions of dynamics, all 
initial conditions and hence all trajectories are considered at the same time, 
forming an ensemble. Two such trajectories are presented on the top panel 
in Fig. ^. For a generic initial point z the trajectory wt goes to infinity 
in the horizontal u direction, while the vertical coordinate vt monotonously 
decreases. However, if the initial point happens to be a point 7(T) on the 
SLE trace, the forward trajectory hits the origin in the mathematical plane 
exactly at time T. 

Conversely, we can fix a point w in the mathematical plane and follow the 
motion of its image Zt under the map ft{w) —£,{t) in the physical plane, with 
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Figure 3: Top panel: two trajectories in the forward flow. Bottom panel: 
the corresponding trajectories in the backward flow. Curved arrows indicate 
the flow of time. In both cases the flow trajectories are shown in grey, while 
the black line represents an SLE trace. The trajectories on the bottom 
panel have been produced with the same noise realization as used for the 
forward evolution, reversed in time as described in the main text. The grey 
trajectories in the top and bottom panel are thus identical. 

the initial condition zq = w. In components zt = xt + iyt, the trajectories of 
this backward flow satisfy the system of equations 

dtxt = 2~~^ ~ dt£,{t), xo = u, 

xf + yf 

dtVt = 2^,* 2 ' = V. (18) 

The two trajectories shown on the bottom panel in Fig. ^ precisely retrace 
the trajectories of the forward flow shown on the top panel. This has been 
achieved by driving the backward evolution ( [l8|) by the time reversed noise 

,^(r — t) — (,{T) = ^(t) as compared to the forward evolution. In this case 
the flnal point zt of the trajectory that started at the origin coincides with 
the tip of the trace 7(r) at that time, but the rest of the trajectory does 
not follow the SLE trace. If we drive the backward flow by an independent 
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copy of £,{t), then even the final point zt wih be different from 7(r), but in 
the statistical ensemble zt and 7(T) will have the same distribution. 

Now we can introduce the probability distribution function of the process 
Zt = xt + iyt in the physical plane defined by: 

p{x,y,t) = {5{xt-x)5{yt-y)). (19) 

From Eqs. ^) it follows immediately that y, t) satisfies the following 
(generalized) Fokker-Planck equation: 



r\, 2x 2y 

dtp{x,y,t) = -dl + c\d^\°' + — 2 ~ ^v^^ — 2 /'(^'y'*)- (2^) 
2 x'^ + x'^ + 



Here (sometimes also written as (— A)"/^) is the Riesz fractional deriva- 
tive, which is a singular integral operator whose action is easiest to describe 
in the Fourier space: if f{k) is the Fourier transform of a function /(x), 
then the Fourier transform of \dx\'^ f{x) is \k\°' f[k). 

As we have discussed, at long times the growth is dominated by the 
stable process in the driving function (^), and we can set k = 0. So our 
main analytical tool is the following Fokker-Planck equation: 

2x ^ 2y 

■ rj.2 y2 y J.2 _^ yZ 



dtp{x,y,t)= c\dX + dx 2 , 2 ~ 2 1 2 (21) 



Let us discuss the boundary and initial conditions for this equation. 
The initial condition for the Fokker-Planck equation ( |2l| ) depends on the 
initial conditions xq = u, yo = v in the stochastic equations (0). For 
the distribution of the SLE tip 7(4) the appropriate initial conditions are 
xq = 0, yo = e, where e is an infinitesimal positive number. For the exact 
Fokker-Planck equation (|20|) this translates into the initial condition 



pix,y,0)=5{x)5{y-e). (22) 



However, for the approximate equation ( piD the situation is more subtle. 
The crossover time to = 0(1)- For t < Iq the drift in the x direction 
(towards x = 0) dominates over the Levy term. For t > Iq the opposite is 
true. A simple picture is then that before the initial 5 function is advected 
by the drift velocity 2/y in the y direction. By the time to it becomes 

po{x,y) = p{x,y,to) = 5{x)6{y-yo), y^ = 2tl'^ ^ c~^'^^-''\ (23) 

This is the initial value that we shall assume for our problem. In the fol- 
lowing sections we will mostly use the notation pQ{x,y), using the explicit 
expression when necessary. 
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Let us comment that if we tried to be more careful and included the 
effects of the Brownian forcing before the crossover time to, then the distri- 
bution at time to would not only be advected to yo but would also broaden 
to a Gaussian with variance kIq. This refinement would not change any 
arguments in the later sections, since all we need there is that the Fourier 
transform in x of the initial distribution is broader than e~'^^^°'^ for long 
times, see the discussion preceding Eq. (|36|). This is a good approximation 
for both the initial distribution (^) or its Gaussian variant for sufficiently 
long times, and becomes better and better as time increases. 

As for the boundary conditions at y = 0, we have no need to be very 
explicit about them, since p{x,y,0) vanishes for y < yo, and our equations 
of motion (|l8|) represent a situation in which yt continually increases as t 
increases, so that p{x, y, t) will also vanish for y < yo at all times t > 0. 

4 Qualitative description, distance scales 

In this section we analyze in qualitative terms the long-time limit of the 
evolution of the tip 7(t), by looking at the consequences of equations ([l^ . 
According to the discussion in the previous section, Re7(t) and lm^{t) have 
the same joint distribution as xt and yt- 

For small times, up to the crossover time to, the drift term in the 
Langevin equation dominates over the Levy noise. Therefore, both xt and 
yt, and ^{t) as well, grow as \/t. For larger times, t ^ Iq, there are two dif- 
ferent characteristic length scales, X{t) and Y{t). In this regime the forcing 
(,{t) is dominated by the Levy process La{t). The probability for a total 
motion X(t) over a time t for this process is described by Eq. (Q). Typically 
the motion is dominated by a single long jump, and the jump has an order 
of magnitude 

X{t) ~ (ct)^/". (24) 

(This can be understood as rescaled fractional moments {\La{t)\^)^^^ , see Eq. 
(^)) Since the typical jumps of ^(t) become arbitrarily large at long times, 
Xt also becomes large, and therefore, the drift term in the first equation of 
( ITsI) becomes negligible. In this limit, xt behaves like the driving force, and 
we find 

t^oo : |xt| ~ A:(t) ~ (ct)^/". (25) 

The Loewner evolution with Levy flights produces, in general, a forest of 
(sparse or dense) branching trees, growing form the real axis. The above 
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Figure 4: Growth of the Levy-SLE parallel to the boundary. Here we plot 
^|^|3a/4^ for various values of a (Brownian motion is set to zero, k = 0). 
The average follows the predicted behavior t^/^. Data collected by averag- 
ing realizations of equation (|l^) for Levy distributed forcing c = 1, time step 
r = 10"^; 10000 runs for a = 0.5, 1.0, 1.3, 3000 runs for a = 0.7, 0.9, 1.1, 1.5. 
The black line is a guide to the eye with the desired slope. The irregular 
points observed in some of the curves are due to large jumps in the forc- 
ing of individuals runs. Such behavior is expected due to the power law 
distribution of the jumps in Levy processes. 



13 



relation then tells us how the forest spreads along the real axis with time. 
This distance is marked out on the plots of trees shown in Figure ^. 

Numerical implementation of the Langevin equation (|l6[), details of which 
are presented in Appendix [A.l|, confirms these qualitative arguments. Fig- 



ure H compares the estimate of Eq. (25) with numerical calculations of the 



trace via simulations of Eq. (16). The agreement is excellent. 

Next we turn to a typical distance Y{t) in the y coordinate. Figure ^ 
clearly shows that this characteristic distance is much smaller than X{t). 
We understand this as follows. If xt were zero, the second equation in (|T8| ) 
would give yt ~ t^^^. Clearly, any non-zero xt only slows down the growth of 
yt- We then conclude that yt, and therefore the height of the trees produced 
by the SLE process cannot grow with time faster than t^/^. Since a < 2, 
it means that Im7(t) always grows slower than Re7(t) and they become 
widely separated at long times. Our major result is that the growing trees 
spread faster horizontally than they grow vertically. Hence, we have 

Y{t) <€. X{t). (26) 

An estimate of the scaling of Y{t) can be obtained from the second 
equation in ([l^) where we replace xt by the Levy process and average over 
it using the probability distribution (^). This gives a typical behavior of yt'- 

dtyt^ rdx^^P{x,t)= rdke~y^\^\~'\'^\"'. (27) 

To estimate the k integral we can drop the term yt\h\ in the exponent, since 
this quantity is of order Y[t)/ X[t) <^ 1. Thus we get 

9m - (28) 

The time integration then gives a result that the length scale for the y 
direction is 

yit) = yo + ^^^^t^--. (29) 

a 

Here yo is formally the constant of integration, but it should really be 
thought of as an adjustable constant inserted to make up for any errors 
we might have made in doing the integrals. In particular, it takes care of 
any effects from the early-time region, where we surely do not have the 
calculation under control. 
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Figure 5: Growth of the Levy-SLE perpendicular to the boundary. We plot 
(y) for various values of a. Same details as in previous figure. Initially, the 
trace grows as \/t for all values of a. This behavior changes around the 
characteristic time to ~ 1- The height of the trace saturates for a < 1, while 
it grows indefinitely for a > 1. This change of behavior demonstrates the 
global implications of the phase transition at a = 1. 
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The phase transition at a = 1 is manifested by a quahtative difference 



between a > I and a < 1. From (29) we can see that while for a ^ 1 the 



average height of the trees grows to infinity as while for a < 1 it 

saturates at a finite value yoo- Figure |5| provides an illustration of the phase 
transition at a = 1 separating different behaviors. More detailed comparison 
between our analytical predictions and numerical simulations is provided in 
the next Section. 

5 Solving the FPE 

In order to quantify these predictions we need to return our attention to the 
Fokker-Plank equation (|2l|). If we perform the Fourier transform in x and 
integrate in time we can find a compact form of this equation, which reads 



pik, y, t) = e-^l^l"*po(A:, y)-dyj dk' dt' e~y\^'\-'\^\''^'~'y{k - k', y, t') 
+ k fdk' f dt' sgn(A;')e-^l*^'l-^l'^l"(*-*')p(fc - k', y, t'), (30) 



where po{k,y) is the Fourier transform of the initial distribution (23). At 
long times p{x,y,t) is spread over the scale X(t) as a function of x. Its 
Fourier transform p{k,y,t), as a function of k, is significantly non-zero on 
the scale X(t)~^. At the same time, due to the exponential factors e"^^^ 
the relevant values of k' in the integrals in Eq. (|30|) are of the order y"^ > 
Y{t)~^. The scale Y{t)~^ is much larger than the range X(t)~^ where p is 
non-zero, hence, when integrating over k' we can use the approximation 



p{k - k',y,t) w 6{k - k') j dk"p{k" ,y,t) = 2Trp{0, y,t)6{k - k'). (31) 
The Fokker-Planck equation then reads: 

KKy,t) = e-'\^\"%{k,y) - 27r / W^l^l-'^l^l"(*-*')a,p(0, y, t') 



JO 

t 



+ A7r\k\ / dt'e-^l^-|-^l'=l"(*-*'V(0,2/,t')- (32) 







This is the main approximation that we will use in order to study the be- 
havior of the Levy-SLE process at large times. 

Notice here that the distribution function p, for every x and t, depends 
only on the initial condition po a-nd the history of the distribution at x = 
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for earlier times t' < t. Therefore, in order to study the probabihty density 
function described by the Fokker-Planck equation, we first need to calculate 
the behavior of this distribution for small x, that is p{0,y,t). Then, by 
substituting in Eq. (|32|), we can in principle estimate the full distribution. 
However, in this paper we are mostly interested in the way this process grows 
in the y direction. Hence, we will first find p{0, y, t) which characterizes the 
growth near x = 0, and then obtain the distribution 



p{y,t)= dx p{x,y,t) = p{0,y,t) (33) 

J — oo 

of y's integrated over all x by setting /c = in (^): 

P{y, t) = po{y) -2Tr [ dt' dyp{0, y, t'). (34) 







This equation immediately leads to the average (y), which is understood as 
the average over all x: 



(y) =yo + 2vr / dt' / dyp{0,y,t'). (35) 
Jo Jo 

Therefore, the distribution and its mean in Eqs. (^) and (^) depend only 
on the behavior at x = at times t' < t. This is a direct implication of Eq. 
(32) and our main approximation (^). 

Let us emphasize again that our approximation works in the long time 
limit. We will assume that we can use approximate expressions in time inte- 
grals for all t > to- Thus, we will treat all time integrals Jq as J^* + correction. 
The corrections come from short times, and we cannot extract them from 
our analysis. They all will be hidden in the terms dependent on the lower 
limit to of the time integrals. In several cases the lower cut-off at tQ is 
necessary to avoid spurious divergencies. 

Let us now consider p{0, y,t). A closed equation for this quantity results 
from integrating Eq. (|3^) over k. To do this we observe that in the first term 
(the initial value at t = 0) for the relevant values of k the function pQ{k,y) 
is much broader in k than e"'^!'''"* at long times. Hence, in the integral over 
k we can replace po{k,y) by its value at /c = 0. Then it follows that 



P(0,.,t) = gM- tdt'piO,y,t')dy- 



2nXit) Jo X{t-t',y) Jo ' 'X{t-t',yy 

(36) 
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where the scale X{t,y) is defined as 

_L_. (37) 

and 

X{t) = X{t, y = 0) = —^^^tV". (38) 

Equation ( |3^ ) is easily solved after performing the Laplace transforma- 
tion in time t. For the transform 

POO 

p(0,2/,A) = / dte-^'piO,y,t) (39) 
Jo 

we obtain an ordinary differential equation 

a„p(„.„ A) + i^lg^KO, A) = 5^?«(0. (40) 



r-oo -At /-oo p-!/Ifc| 

K{X,y)= I dt——= dk-——, (41) 



where 



X{t,y) A + c|A;|°' 

and -fC(A) = K{X,0). Using the initial condition po(0,y) = S{y — yo), the 
straightforward solution of Eq. (^) is 

,„ J. K(\) K{X,y„) I f'_dy_\ „„, 

The inverse Laplace transform of this solution gives p{0,y,t). 

Notice that ( ^2| ) is valid only for y > yQ. Since our approximations 
only work at long times, we expect our solution to give good results for 
y ^ yo- The approximations will usually result in the necessity to introduce 
a fitting parameter (called "correction" in the discussion after Eq. (]35|)) in 
the time evolution of averages for the process. Moreover, there is an upper 
cut-off that stems from the Langevin equation and the fact that y cannot 
grow faster than t^^'^ (see previous section). Since we used this fact while 
making the approximations that lead to Eq. (^), the range of validity of 
our solution is yo ^ y ^ 

In the following we will analyze the properties of the distributions p(0, y, t) 
and p{y,t) in three separate cases a > 1, a = 1 and a < 1. For each case 
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we will repeat the following steps: first we calculate p{0, y, t) from Eq. 
then, by substituting this solution into Eq. (|3^), we will calculate the av- 
erage height (y) and the distribution p{y,t). In these calculations we need 
approximate expressions for the function K{X,y). These expressions are 
derived in Appendix I 



5.1 Results for a > 1 



In this case we can use the approximation (|8q ) from Appendix A^ for K{X) 
and K{X,y). Eq. (^) then gives 

1 ' iAi-l/<.,-> ■ 2^ 



p(0,„,A)»-exp^--V-"^j, A = -^j^. (43) 

a 

To calculate the time dependence of the distribution we take the inverse 
Laplace transform: 



a+oo A \ 

2m ' ^ ' 



As usual, the integration contour in the last equation goes along a vertical 
line Re X = a, where a should be greater than the real part of any singularity 
of the integrand. Changing the integration variable to At we obtain that 
answer which, apart from the overall prefactor 1/t, has acquired the form 
of a scaling function: 

p(0...,)»^F(S), s.^. r«) = J^^,-i (45) 

Since the scaling function F{y) depends only on the combination yt~^'^^^°', 
its derivatives with respect to y and t are related: 

(T f 

dyF{y) = --dtF{y). (47) 

a — ly 

The integrand in Eq. (^) contains a branch cut which we choose to run 
along the negative real axis. The integration contour can be deformed to go 
from — oo to along the lower side of the cut, and then from to — oo along 
the upper side. This leads to the final answer for the scaling function F{y): 

F{y) = i TdAe-^-l-^^l^^^^'-^sinf sin-Ai-i/"yV (48) 
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y/Y(t) 

Figure 6: The distribution of heights scales as y/Y{t), where Y{t) is given 
by Eq. ( ^5[ ) for a = 1.3. The distribution is shown at three different times 
(black, red and green curves), all within the limiting region of large times 
where asymptotic behavior y oc holds. 



The overall prefactor 1/t in p{0,y,t) can be understood as follows. The 
distribution p{x,y,t) at long times spreads in the x direction up to scale 
X{t), and in the y direction up to scale Y(t). The total area "covered" by the 
distribution scales with time as X{t)Y{t) oc t. Therefore, at the particular 
value X = the density p(0, y, t) decays with time as 1/t. However, if we are 
looking at the distribution of the y coordinate for x = 0, and its moments 
(?/"■), we should divide p{0,y,t) by the normalization 







dypiO,y,t) = (49) 



m 

2^tr(l-i) 

The normalized distribution is then 



Pn{0,y,t)^^^^^y^F{y). (50) 

Moreover, the integrated distribution p(y, t) exhibits the same scaling as 
p{0,y,t). Indeed, using the relation (^) in Eq. (U) we obtain: 

1 a I 

piy,t)=Poiy) + T^ T-F{y). (51) 

Y{t) a-ly 
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y/Y(t) 



Figure 7: The distribution p{y,t) as a function of y/Y{t). The theoretical 
prediction Eq. (^) (soUd curve) and the numerical distribution (black dots) 
are different, however, they have a similar dependence on y for the values 
where we believe the solution is valid, t^"^ <^ y <^ t^l'^ . Here, Y{i) = 50, 
tV2 ~ 300 and t^^ 1, so the region of validity of Eq. (|5^) in the scaled 
variable is 0.02 <C y/Y{t) <C 6. This explains the disagreement between the 
theory and the numerics for y/Y{t) > 2. Also, for small values of y/Y{t) 
where we should not trust Eq. (|5l|), the theory still gives a significant 
weight to the distribution p{y,t). This is, presumably, the reason for the 
discrepancy between the numerics and the theory in the range y/Y{t) < 2. 



Fig. ^ shows the scaling collapse of the numerically calculated distributions 
p{y,t) for a = 1.3 and three different times. We see that, indeed, p{y,t) is 
a scaling function of y/Y[t) in agreement with our predictions. 

We can calculate the asymptotics of the function F{y). For small values 
of y we can neglect the term with y in the exponential in Eq. (|4^), as well 
as replace the sine function under the integral by its (small) argument: 

ct — 1 1 , . 

F{y « 1) « —j-y. (52) 



a 



For large y we need to use the steepest descent method for the contour 
integral in Eq. (^), which results in 



Fiy » 1) 



(3 



1/2 



a 



a 



1 \"/2 

-y 



exp 



a 



1 



a 



1 



a 



-y 



(53) 
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Figure 8: The average height {y) for SLE driven by Levy flights with a = 1.3 
grows as a power-law The red dashed line is a fit to Eq. ( |54|) for 

t > 1, where we only vary the parameter yg- 



We have to remember that we can only trust this result for yo y t^l"^ . 

Figure |^ shows a comparison between the numerical data and the theo- 
retical prediction of Eq. ( |5l|) for the distribution p{y,t). While the overall 
dependence on y is similar between the two, we would obtain a better fit 
for yo ^ y ^ t^^"^ if we redistributed the weight outside this region to the 
range were Eq. (51) is valid. 

Next, we will calculate the time evolution of the average height of the 
growing trees {y) from Eqs. (35, E^): 



{y) = yo + 



r(i 



yo + 



2 ri 



a/a I 



(54) 



Here, all short time contributions are included in yo- This nicely fits the 
numerics, see Fig. |8|, and reproduces the result (29) of the simple argument 
using the Langevin equation. 

We also want to compare the distribution at x = to the distribution 
averaged over all x. We calculate the average value (y)° (the superscript 
indicates that this average is calculated at x = 0) from the distribution 

r 



ac 



l/a 



IT 




In 




COS — 


r(:- 






a 




a/ 





(55) 
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0.001 0.01 0.1 1 10 

x/X(t) 

Figure 9: Average height of the trace y = Im7(t) as a function of x/^ct)^^" 
for SLE driven by Levy flights with a = 1.3. y data are bined logarithmically. 
Here, (ct)^/" = 1827.15, 6217.82, 21156.6 for the three values of time. The 
average height close to x = is 2 times bigger than the height at large x 
and roughly 1.4 times higher than the global average {y). The theoretically 
predicted value of the ratio between the height at x = and the average is 
1.9 (Eq. (|5^)). This discrepancy is most probably due to a finite time effect 
and the limited amount of data close to a; = 0. 
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The ratio of the two averages (neglecting yg) is 



27r 
sm — 

a 



This tends to 2 as q — > 1 from above, and to 7r/2 as q — > 2 from below. We 
observe similar behavior in our numerical results, where the average of y at 
x = is higher than the overall average (Fig. However, the ratio ( |56| ) 
is not matched exactly. Presumably, this is because we do not have enough 
data close to x = and we cannot reach long enough times in order for the 
various constants (like i/o) to be negligible, so that Eq. (^) is accurate. 

5.2 Results for a = \ 

Now we use the approximations (|7^, from Appendix A. 2 in Eq. 



The resulting expression for /9(0, y. A) is difficult to analyze without further 
approximations. We will evaluate it as well as its inverse Laplace transform 
with logarithmic accuracy, which amounts to three assumptions. First, we 
assume that all the logarithms that appear are large compared to constants 
of order one such as vr, c, etc, which will be neglected. Secondly, the loga- 
rithms are assumed to be small compared to power laws for large arguments: 
Ini ^ t. Finally, the logarithms are slow functions as compared to power 
laws and exponentials, and in integrals can be replaced by their values at the 
typical scale of variation of the fastest function under the integral. All sub- 
sequent equations in this section will be obtained with logarithmic accuracy 
using these assumptions. 
First we have 

„„,,,A).±;<|Se.p(_£_|). ,57) 

The time dependence now follows from the inverse Laplace transform, using 
the same contour integral described in the previous section: 



p(0,y,t)pa J- rd\ ^ Ato^'' A^o ( ^_^y_\ exp f- A - - ^ 
^ ^ In^f V21n2|t; 2 In || 



(58) 



The integral of this expression over y 

. ,^ . 1 'L ^ 



/ dyp{Q,y,t)^— (59) 
Jo TTci In^ (%^) 
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time 



Figure 10: The average height {y) for SLE driven by Levy flights a = 1.0 
grows logarithmically with time. The red dashed line is a one parameter fit 
for t > 10 to the predicted function A + | Int. 



leads to a normalized distribution at x = 0: 
c lnM4^) 



2-K In f In ^ 

to yo 



. [tt cy \ ( . c y 
dX ^ — -r-^ sm ^ — -r exp — A r 



(60) 



The mean value of the height of trees near x = follows from yOn(0, y, t) 
using the same arguments as before: 

(y)0«^ln^. (61) 



The average height (over all x) is also found easily from Eq. (|35[): 

2 /■* (If' In 1^ In — 2 



As shown in Fig. 10 this is in good agreement with the numerics. The ratio 
of the two averages in the long time limit is {y)^ / {y) = 2, consistent with 
the hmit a ^ 1 of Eq. (|5|). 
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The asymptotics of the distribution pn(0, y, t) for small and large values 
of y/lnct can be found similar to the case a > 1: 



fc2 In^J^ 

4 ct 



Pn{0,y,t) « < 



y 



4 In^ 

3/2 1„2 In — 

c/ In — m^^^^m^^^ 
47rV2 In ^ In ^ In^ ^ ^ ^ 



to yo 



cy 



21n^ 



, Inct < y < t^/^. 
(63) 



Finally, using Eq. (|34|) , we get an expression for the integrated distribu- 
tion: 



p{y, t) w po(y) + — In — / dA 
2vr to Jo 



t In ^ In ^ 



J 3 



sm — 



vr cy 



At/ 



exp — A — - 



c y 



21nf^ 
(64) 



The asymptotics of these expression follow as before: 

f In^ 

to yo 



fcSln^f In-^* 



Piy,t) -po{y) ~ < 



4 In 



5 ct 



-y, 



8t i„ 8t 



c3/2 In fin -2^ In 



47rV2 



y^' exp 



cy 



21n8* ^' 



- ^ 



yo < y < inct, 

Inct < y < t^/^. 
(65) 



5.3 Results for a < 1 

In this case we use the approximation (|8^) from Appendix A. 2 leading to 

,2-a 



p(0,y,A) «iVK(A)y2-2"exp( 



C=-r(l-a), N 

c 



y 



C{2 - a) 
27rC' 



(66) 
(67) 



The inverse Laplace transform of this expression gives the leading approxi- 
mation 



p(0,y,t) 



N 

W) 



y2-2-exp 



(l-")c_^2-a 



2r(3 - a)' 



(68) 



The obtained result depends on time only through the overall factor 
X~^{t). We can understand this as follows. The distribution p{x,y,t) at 
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Figure 11: Distribution fits for a < 1. We compare the numerically calcu- 
lated distribution p{y, t) to the theoretical curve for p^o given by Eq. (|69|) 
with one free parameter for normalization. We claim that p{y, t) = /9oo(0, y) 
for yo <^ y <ti t^^^ where our solution is valid. 



long times spreads in the x direction up to the scale X{t) but becomes 
stationary in the y direction. Therefore, at the particular value x = the 
density p{0,y,t) decays with time as X~^{t). However, if we are looking at 
the distribution of the y coordinate for x = 0, we should normalize Eq. ( |68| ) 
which gives the truly stationary distribution (normalized by the appropriate 
choice of A^i) 

.2-2a / (1 - a)c „,2- 



pUO, y) « A^iy^-^" exp ( - ^^^^J^f '" ] , (69) 



in agreement with numerics, see Fig. 11, where we actually observe that the 



integrated distribution coincides with Poo{0,y) at long times. 

The stationary distribution (|69| ) allows us to calculate the average satu- 
rated height of the trees: 

yoo = J^ dyypoo{0,y) = [ ) r-i(^-— jc-i/^^-). (70) 

This is in very good agreement with the numerically calculated values shown 
in Fig. IT2. 
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Figure 12: (Left) The ratio of the numerically calculated over the theoret- 
ically predicted value of the saturated height Uoo for a < 1. (Right) The 
saturated height of the trees vs. the strength. The dashed lines are the 
theoretical values for yoo, Eq. (|70|). 




Figure 13: Distribution of heights averaged over all x for SLE driven only 
by Levy flights and a = 0.7. The distribution is shown at three different 
times (black, red and green curves), corresponding at the limit of large time. 
We see how the distribution of the height of trees is stationary. 
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Figure 14: The average height (y) for SLE driven by Levy flights a = 0.7 
saturates to {y)oo as t^~^/°' . The red dashed line is the analytic result, Eq. 
([7l|), with the value of = J^i^jt^^J^^ obtained in Eq. (29). (y)oo was 
calculated from the two numerical points of y for the largest times. 



Let us discuss now the integrated distribution p{y,t) and its mean (y). 
Unfortunately, in the present case {a < 1), the Eqs. ( p^ ) and (^) do not 
give reliable results simply because the apparent distribution and saturation 
height are very sensitive to the lower limit to, and the results are of the 
same order as the initial conditions at to. Analytically, we can see that the 
distribution p{y, t) becomes stationary as t — > oo, even though we cannot 
determine oo). The time independence of the distribution p{y,t) at 
long times is checked numerically in Fig. |l^. Numerics presented in Fig. 
11 indicate that p{y,t) = poo(0,y) (see Eq. (pgj)) for the appropriate range 
yo <^ y <^ t^^'^, and we will discuss why this is true below. 

We can also see that the way the average tree height approaches its 
limiting value is given by the power law 



(y) = (y)oo-I?c-V"ti-i/". 



(71) 



We have previously calculated D in Eq. ( p9| ) using the Langevin formulation 
of the process. This result agrees well with numerics, as demonstrated in 



Fig. |14. 

We can argue that 



yoo and that p{y, oo) 



Pc 



if we return to 



the initial description of the process seen as SLE trees growing forward in 
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Figure 15: Average height of the trace y = Im7(t) as a function of 
x/{ctY/°' for SLE driven by Levy flights with a = 0.7. y data are 
binned logarithmically and the average of every bin is plotted, (ct)^/" = 
2.9 10^2.8 10^,2.75 10^. 



time fT]. For a < 1 the jumps of the Levy process are large and we know 
that any new tree is most likely to grow starting from the real axis. The 
trees are sparse and the new tree will grow isolated from its neighbors, hence 
it will be identical to any other tree, including the trees that grow close to 
the origin at x = 0. Therefore, we expect the distribution of y's at x = 
to be identical to the distribution at any other x. Numerics also support 
this argument. In Fig. O we observe that the average height of the trees 



is practically independent of the value of x. Also, in Fig. [llj we compare 
p(y, oo) and poo while in Fig. 12 we show that {y)oo = Voo- 



6 Conclusions 

In this paper we have analyzed the global properties of growth in the complex 
plain described by a generalized stochastic Loewner evolution driven by a 
symmetric stable Levy process La{t), introduced in our previous paper [j^. 
The phase transition at a = 1 whose implications for local properties of 
growth were the subject of Ref. Q, also manifests itself on the whole plane 
resulting in a rich scaling behavior. 
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We have used a Fokker-Planck equation to study the joint distribution 
p{x, y, t) for the real and imaginary parts of the tip of the growing trace. 
The presence of the Levy flights in the driving force imposes very different 
dynamics in the x and y directions. While in the x direction the process 
spreads similarly to the Levy forcing x ~ X{t) ~ t^/", the SLE dictates 
y <^ X{t), for all values of a. This separation of the horizontal and vertical 
scales in the process allows us to make sensible approximations and explore 
geometric properties of the stochastic growth in all phases, a < 1, a = 1, 
and a > 1, both qualitatively and quantitatively. 

For a < 1, the vertical growth saturates at a finite height y^Q. In terms 
of the picture presented in Q, long jumps occur often so that new trees 
grow isolated and there is a small chance that the trace grows on an already 
existing tree. 

For a > 1, the average height of the process grows as a power law t^-i/" 
with time. New trees grow close to old ones, so that when the process 
returns to a previously visited part of the real axis it will have to grow on 
top of already existing trees. Eventually the trace will grow past any point 
on the plane. 

At the boundary between the two phases, a = 1, the height of the process 
grows logarithmically with time. 
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A Appendices 

A.l Numerical calculations 

The interpretation of equation ( [l6| ) is very helpful to our calculations. Zf 
and the tip of the trace have the same distribution. This allows, instead 
of calculating the trace 7(t) for every time t and noise realization (0(n^)), 
to efficiently collect statistics for the position of the tip by integrating the 
Langevin equation (|l^) (0(n)). 

Following Ref. we approximate ^{t) by a piecewise constant function 
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with jumps appropriately distributed: ^(t) = for (j — 1)t < t < jr. For 
such a driving function the process zt in Eq. ( p!6|) can then be calculated 
numerically as an iteration process of infinitesimal maps starting from 
the condition z = as follows: 

Zn = z{nT) = /„ O /„„! ... o /i(0) - in- (72) 

The infinitesimal conformal map /„ at each time interval n is defined by: 

fn{z) = W-\z) = ^{Z - " 4r + (73) 

The value of f„ is randomly drawn from the appropriate distribution. The 
number of steps necessary to produce an SLE trace up to step n grows only 
as 0{n). All numerical results in the next section have been calculated using 
the average of Eq. ([T^ ) over many noise realizations. 

The trace can also be produced directly IQ, as g~^{i{t),t), in which case 
we approximate 

7j = jUt) = /i . . . o o fniU)- (74) 

However, the number of steps in this method grows as O(n^). We used 
this method to verify that numerically calculated z and 7 have identical 
distributions. Eq. (|7^ was also used to calculate the traces shown in Fig. 

Here, we will assume k = for simplicity, that is, the driving force is 
pure Levy flights .^(t) = c^^°'La{t). The addition of a Brownian motion will 
not affect our conclusions. For all realizations of the Levy-SLE process we 
take c = 1 and r = 10~^ unless otherwise noted. 



A. 2 Asymptotics for K{\^ y) 



Let us consider (we need to use the lower cut off to here to have a convergent 
result for a < 1) 



K{\) 



00 g~-^* 



dt 



to 



Xit) 



2r(i + ^: 



dtr^^^e 



xt 



to 



( 2Vil + ^"l 



^ c 



a = 1, 



(75) 



where r(a, x) is the incomplete gamma function, and Ei{x) is the exponen- 
tial integral. Since A has the dimension and the meaning of frequency, and 
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we are interested in t ^ to, we will only need the small argument asymp- 
totics of these functions: 

T(a,x)Kir(a) , Ei(x)^-lnx, x <Cl, (76) 

a 

This gives for Mq <^ 1 

K{X) « ^A-i+^/" + Btl'^^", a / 1, (77) 

A=—^ , B = 4} ^r(l + i), (78) 

acV" sin ^' cV" l _ « ^ ^ ^ 

a 

K{X)^-ln^, a = l, (79) 

C AlQ 

For Q > 1 we can set tQ = and obtain 

K{X) = AX-^+^Z'', a>l (80) 

and for a < 1 we can set A = 0: 

K{0) = Btl~^^", a<l. (81) 

We now turn to the Laplace transform K{X,y): 

'^(^•'') = Vo "'XT*^- '''' 

Since in the Laplace transform the important values of A are the inverse 
typical time scales, this means that the relevant asympotics of K{X, y) are 
those with Xy°'/c <^ 1. The opposite case of Xy'^/c ^ 1 corresponds to 
short times, where our basic approximation is invalid. So from now on we 
will focus on the limit Xy°'/c <C 1. 

This integral can be evaluated exactly in a number of cases. First, when 
y = 0, the integral converges for a > 1 and gives the same expression as 
K{X) in Eq. (|80|). Secondly, for A = the integral converges (for y > 0) for 
a < 1 and gives then 



2 2 
K(0, y) = - dk k-'^e-y^ = Cy'^'^ , C = -T(l 
c Jo c 



a). (83) 



All the constants A, B, and C defined above diverge as l/(a — 1) as a — > 1. 
Finally, for a = 1 we get 

KiX,y) = -e'y/^E,{^) « ^Inf , ^ « L (84) 

c c Xy c 
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In general for Xy'^ /c ^ 1, a good approximation for K{\,y) is the sum 
of expressions in Eqs. (180, |83|): 



i^(A, y) « AA-i+i/" + Cy"~^ (85) 



Not only this approximation reproduces the correct limits in Eqs. ( |80D and 
( p3D , but in the limit a ^ 1 it also reduces to Eq. ( p^ ) . This approximation 
can be obtained by splitting the k interval in the integral in Eq. ( ^2|) into 
two at the value ko = (A/c)^/" and in each resulting integral replace the 
denominator by the largest term in it. 

Notice that for a > 1, and in the limit of interest Ay°/c <C 1 the first 
term in Eq. ( ^5|) is much greater than the second, and we can use Eq. (|80| ) 
for both K{X) and K{X,y): 

K{X) « K{X, y) « ^A-^+^/°, a > 1. (86) 

For a < 1 the opposite is true, and we can use Eq. (83) as a valid approxi- 
mation: 

K(A,y)«Cy""S C = ^r(l-a). (87) 
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